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Delay  transition  using  non¬ 
equilibrium  C02 


CALtit 


PROBLEM:  In  hypersonic  flight,  heating  loads  are  typically  a  dominant  design  factor 

Turbulent  heat  transfer  rates  can  be  about  an  order  of  magnitude  higher  than  laminar  rates  at 

hypersonic  Mach  numbers 

A  reduction  in  heating  loads  by  keeping  the  boundary  layer  laminar  longer  means  less  thermal 
protection  needed  and  hence  less  weight  to  carry,  or  conversely  more  payload  deliverable  for  a 

given  thrust. 

OBJECTIVE:  Delay  transition  from  laminar  to  turbulent  flow  in  the  boundary 
layer  of  a  slender  hypersonic  body  by  using  nonequilibrium  C02 


Transition  in  high  Mach  numbers 
occurs  through  the  Mack  mode  - 
amplification  of  acoustic  waves 
traveling  in  the  boundary  layer 


Molecular  vibration  and 
dissociation  damp 
acoustic  waves 


^ 


At  relevant  conditions,  C02 
absorbs  most  energy  at  the 
frequencies  most  strongly 
amplified  by  2nd  (Mack)  mode 


Inject  C02  to  delay  transition  in  air  flows  of  interest 
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Background 


GALffl 


Experimental  data  show  that  transition  is  delayed  for  C02  flows  compared  with  N2  and  air 
flows  for  a  given  stagnation  enthalpy,  h0 

These  observations  point  to  a  second  mode  transition  (or  Mack  mode)  for  the  conditions 
studied  as  well  as  to  the  importance  of  non-equilibrium  effects  of  C02  on  stabilizing  the  flow 


0  5  10  15 


h0  (MJ/kg) 


C02  Transition  Re*  is 
about  5X  that  of  Air  and  N2 


From  Hornung,  H.G.,  Adam,  P.H.,  Germain,  P.,  Fujii,  K.,  Rasheed,  A.,  “On 
transition  and  transition  control  in  hypervelocity  flows,”  Proceedings  of  the  Ninth 
Asian  Congress  of  Fluid  Mechanics,  2002 
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Background 


GALM 


Computations  show  that  when  pure  C02  is  in  vibrational  and  chemical  non-equilibrium, 
these  relaxation  processes  absorb  energy  from  acoustic  disturbances  in  the  boundary  layer 
whose  growth  is  responsible  for  transition  in  hypervelocity  flows 


Confirms  trends  seen  in  experiments  where  C02  exhibits  delayed  transition  with  respect  to 
Air  or  N2  for  h0~5-10MJ/kg 


For  air  -  no  effect  from  vibrational  relaxation 
and  chemical  reactions  on  stabilizing  the 
boundary  layer 


ForC02_  vibrational  relaxation  and  chemical 
reactions  stabilizes  the  boundary  layer 


co2 


From  Johnson,  H.B.,  Seipp,  T.G.,  Candler,  G.V.,  “Numerical  study  of  hypersonic  reacting 
boundary  layer  transition  on  cones,”  Physics  of  Fluids,  10  (10):  2676-2685  Oct.  1998. 
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Amplification  rate/cycle  or  absorption  rate/wavelength 
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Background 


GALffl 


Computed  acoustic  absorption  rates  (open  symbols)  -  Fujii  et.  al 

Computed  acoustic  amplification  rates  (solid  symbols)  -  after  Reshotko/Beckwith  and  Mack 
For  C02  the  broad  sound  absorption  curve  peak  coincides  with  the  amplification  peaks 
This  coincidence  is  most  pronounced  at  enthalpies  of  ~10  MJ/kg 


From:  Fujii,  K.,  Flornung,  FI.G,  "Experinmental  Investigation  of 
High-Enthalpy  Effects  on  Attachment-Line  Boundary  Layer 
Transition,"  AIAA  Journal,  Vol.  41,  No.  7,  July  2003 
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Facility:  T5  Hypervelocity  Shock  Tunnel 


Impulse  Facility,  test  time  in  the  order  of  ms,  but 
high  stagnation  enthalpies  and  pressures 
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Free-stream  mixtures  with  C02  CAL^ft 


From  Leyva,  IA,  Laurence,  S,  Beierholm,  AK-W,  Hornung,  HG,  Wagnild,  R,  and 
Candler,  G,  “Transition  delay  in  hypervelocity  boundary  layers  by  means  of 
C02/acoustic  instability  interactions,”  47th  AIAA  Aerospace  Sciences  Meeting, 
AIAA  2009-1287. 
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Porous  Injector 


GALM 


Porous  Injector  Rationale 


Move  to  a  transpiration-like 
approach  instead  of  discrete 
jets 

High  velocity  jets  disturbed  the 
boundary  layer  -  penetrate  to 
shock  layer 

Need  lower  flow  penetration  into 
the  boundary  layer 
Can  achieve  same  flow  rates  as 
with  jets  -  0-50  g/s 


Porous  injector  section 
Sintered  316LL  Stainless  Steel 
10  pm  media  grade 
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Porous  Injector  Results  (10  MJ/kg) 


GAL&t 


10-micron  Porous  Injector  (Ar  injection  at  3.7  grams/sec) 


Heatflux  distribution  -  Shot  #2600 


-0.9  4.8  -0.7  -0.6  -0.5  -0.4  -0.3 


Cone  surface  coordinate  (m) 


Heat  Flux  (MW/m2) 


Average  over  test  time 
(1.500-2.100  ms) 
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Porous  Injector  Results  (10  MJ/kg)  GAL^ft 


10-micron  Porous  Injector  (no  injection) 


Heatflux  distribution  -  Shot  #2598 


0.5  1  1.5  2  2.5  3  3.5  4  4.5  5 

Heat  Flux  (MW/m2) 


-0.9  -0.8  -0.7  -Q.G  -0.6  -0.4  0  3 


Cone  surface  cccrdinjl*  (ml 


Htai  Fluy 


Average  over  test  time 
(1.500-2.100  ms) 
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Porous  Injector  Results  (10  MJ/kg)  GAL$tt 


10-micron  Porous  Injector  (C02  injection  at  3.7  grams/sec) 


Heatflux  distribution  -  Shot  #2590 


Htaifluit  diamfcuiiMi  -  Sh6T  ff£59G 


0.5  1  1  5  2  2.5  3  3.5  4  4.5  5 


Hut  Flux 


Average  over  test  time 


(1.500-2.100  ms) 
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Transition  Reynolds  Number 


Porous  Injector  Results  1/2  CAL$t 


Re  vs.  injection  <5>  ~10  MJ/kg,  ~55  MPa 


5.50E+O6 

5.00E406 

4.50E4O6 

4.00E406 

3.50E406 

3.00E+06 

2.50E+06 


0  2  4 


41 


6  8 


'  C02 
■  Argon 
Smooth 


mass  flux  (g/s) 
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Stanton  number  (St)  (-] 


Data  reduction:  Average  Heat 
Transfer  Method 


CALtit 


Plot  of  St  vs  Rex  for  T5-2598:  P0  =  55.3  MPa.  h0  =  10.26  MJ/kg,  T0  =  6204  K' 


Plot  of  St  vs  Rex  for  T5-2600:  P0  -  54.7  MPa.  h0  =  9.88  MJ/kg,  T„  =  6044  K 


No  injection:  porous  section 


Plot  of  St  vs  Re*  for  T5-2590:  P0  =  55.9  MPa.  h0  =  9.92  MJ/kg,  T0  =  6064.3  K 

-2 


Argon  Injection:  3.7  g/s 


C02  injection:  3.7  g/s 
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Intermitten cy  (7)  Intermittency  (7) 


Data  Reduction:  Intermittency  Method  G  Lilt 


Shot  2596  Intermittency  (Average  by  Gauge  Position) 

0.25 

0.2 

0.15 

0.1 

0.05 

0 

0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1 


Shot  2600  Intermittency  (Average  by  Gauge  Position) 
1 

0.8 
0.6 
0.4 
0.2 
0 

0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1 


x  x 


X  x  X 


X  XX 


X  x  X 


X  X  X  X  X 


Turbulent  intermittency 


•  Mee  D. J.  and  Goyne  C.P.  (1996)  Turbulent  spots 
in  boundary  layers  in  a  free-piston  shock  tunnel 
flow.  Shock  Waves,  Vol.  6,  No.  6:337-343. 

•  Narasimha  R.  (1985)  The  Laminar-Turbulent 
Transition  Zone  in  the  Boundary  Layer.  Progress  in 
Aerospace  Sciences,  Vol.  22,  29-80. 


Alternate  method  to  determine  transition  location 
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Transition  Reynolds  Number 


yi-*s  ■ 


Porous  Injector  Results  (10  MJ/kg)  GAL&it 


Re  vs.  injection  @  ~10  MJ/kg,  ~55  MPa 


♦  C02 

•  Argon 

>  Smooth 


mass  flux  (g/s) 


Summary  of  results 
(intermittency  method) 
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Theoretical  Injection  CAL&t 


•  Free-stream  gas  is  Nitrogen 

•  Different  injection  geometry 

—  5  degree  cone,  transpiration  from  10  to  90  cm  on  the  cone 


Stagnation  Conditions 

Free-stream  Conditions 

Pressure  (MPa) 

55.0 

Density  (kg/m3) 

0.051855 

Temperature  (K) 

6958 

Temperature  (K) 

925.5 

Entalpy  (MJ/kg) 

9.39 

Velocity  (m/s) 

4039.7 

•  Injection  based  on  profile  suggested  by  Malik* 

—  mdot  based  on  edge  conditions  and  parameter,  fw 


A/2Re  p  v 

_  V  x  r  w  w 

•  For  these  cases,  fw  held  constant  J »■  ”  ~~ 

re  Ve 

—  Mass  flux  decreases  down  the  length  of  the  cone 


*  Malik,  M.  R.,  "Prediction  and  Control  of  Transition  in  Supersonic  and  Hypersonic  Boundary  Layers," 

AIAA  Journal,  vol.  27,  no.  11,  November  1989,  pp.  1487-1493. 
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fw 

Mass  Flux  (g/s) 

0.05 

1.24 

0.1 

2.47 

0.2 

4.96 

0.3 

7.43 

0.4 

9.92 

0.6 

14.88 
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Computational  Results  GALffit 


•  Transition  occurs  at  N  =  ~9.2 

•  Significant  transition  delay 
vs.  smooth  cone  for  C02 

•  Air  and  N2  injection  both 
promote  transition 

•  Mass  flux  =  2.5  grams/sec 
(but  over  entire  surface; 
“Malik  cone”) 

•  Smooth  cone,  xtr  =  63  cm 


N  factor  versus  distance 
along  cone  surface 


From  Wagnild  et  al  2010 
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N  Factor 


Computational  Results  CALM 


Transition  delay  predicted 

—  Increase  in  C02  initially  results  in 
further  delay  but  further  increase 
causes  more  amplification 

For  Ncr  =  9.2 


Pre-heating  further  delays  transition 

—  C02  able  to  absorb  acoustic  energy 
earlier 

—  Higher  temperature  gas  could  also 
contribute 

For  Ncr  =  9.2 

—  Smooth  cone,  xtr  =  63  cm 

—  fw  =  0.1,  xtr  =  83  cm 


Smooth  cone,  xtr  =  63  cm,  fw  =  0.1, 

xtr  —  72  cm  From  Wagni|d  et  g|  2010 
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N  Factor 


Effect  of  gas  and  Temperature  on  Transition .  hlf 

Delay:  CFD  predictions  O/ILCr/ 


0.2  0.4  0.6  0.8  1 

s  (m) 

Alternate  gases  only 
increase  disturbance 


Smooth  Cone 
C02 1*  =  0.1 
Airfw=  0.1 
N?fw  =  0.1 


30 


x 

C  20 


0 

D) 


O 


0 

O 

k_ 

0 

CL 


-10 


-20 

CFD  predicts  that  for  the  current  porous 
design  and  longer  porous  injectors 
transition  could  be  delayed  for 
optimum  flow  rates  and  temperature 
of  C02 


i _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ L 

0  5  10  15 

mdot  (g/s) 


Distribution  A:  Approved  for  public  release;  distribution  unlimited 


19 


Citations 
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•  Aerospace  America  2009  -  Year  in  Review:  Fluid  Mechanics 

“AFRL,  Caltech,  and  the  University  of  Minnesota  have  collaborated  in  a  numerical  and 
experimental  study  on  control  of  high-speed  boundary  layers.  The  team  has  demonstrated 
significant  delays  in  transition” 

•  Annual  Reviews  of  Fluid  Mechanics  2011,  43:79-95.  -  Federov,  A. 
“Transition  and  Stability  of  High-Speed  Boundary  Layers”: 


“..Another  way  to  stabilize  the  second  mode  and  thereby  delay  transition  is  to  add  C02  into  high  enthalpy 
boundary-layer  flow  (Leyva  et  al.  2009).  The  motivation  for  this  new  technique  lies  in  the  following  findings:  Molecular 
vibration  and  dissociation  suppress  the  acoustic  instability,  and  at  relevant  conditions  for  hypersonic  flight,  C02  absorbs 
energy  most  strongly  in  the  frequency  band  associated  with  the  second  mode. 

The  experiments  of  Leyva  et  al.  (2009)  on  a  sharp  slender  cone  in  the  GALCIT  T5  tunnel  showed  that  the  C02/N2  free- 
stream  blends  (without  injection)  lead  to  significant  delay  of  transition.  The  transition  Reynolds  number  more  than  doubled  for 
mixtures  with  40%  C02  mole  fraction  compared  with  the  case  of  100%  N2.  A  similar  effect  was  noted  in  experiments  using 
mixtures  of  air  and  C02  as  the  test  gas.  Experimental  and  numerical  studies  of  the  C02  injection  system  suitable  for  this 
LFC  concept  are  in  progress.  The  effect  of  the  injection  and  the  transition  location  is  gauged  by  solving  the  PSEs  and  using 
the  semiempirical  eN  method  (Wagnild  et  al.  2010)..” 
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Stanton  number  (St)  {-] 


. «!>■  Half-Porous/Half-Smooth  Injector  r.AI  fit 

■  cn*«4rvin  ■  i  emvim  ^  f  H  ■  j  mfj  Mi. 


Plot  of  St  vs  Re*  for  T5-2657:  P0  =  78.3  MPa.  h0  =  8.6  MJ^kg.  T0  =  S577.5  K 


Reynolds  number  (Re^  [-]  p|ot  0|  st  vs  Re,,  for  T5-2653;  P,  =  64.8  MPa.  h0  =  8  MJ/kg,  T0  =  5278.9  K 
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Latest  -> 


Boundary  Layer  Temperatures 


CALtit 


Heatflux  distribution  *  Shot  #2582 


-0.9  -0.8  -0.7  -0.6  -0.5  -0.4  -0.3 

_ Cone  surface  coordinate  (m ) _ 


0.5  1  1.5  2  2.5  3  3.5  4  4.5  5 

Heat  Flux  (MW/m2) 


H0  (MJ/kg) 

P«  (MPa) 

T*  (Eckert) 

T 

Aedge 

5.69-5.92 

~  30 

1465  K 

779  K 

8.29 

53.8 

2149  K 

1272  K 

8.60 

78.3 

2246  K 

1355  K 

8.85 

53.2 

2295  K 

1395  K 

9.46-10.32 

~55 

2725  K 

1728  K 
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Stanton  number  (St)  [-) 


Ongoing  Challenges  GALm 


Plot  of  St  vs  Rex  for  T5-2598;  P0  =  55.3  MPa,  h0  =  10.26  MJ/kg,  T0  =  6204  I 


10.3  MJ/kg,  55  MPa 
Natural  transition  @  72  cm 
Delay  observed 


8.6  MJ/kg,  78  MPa 
Natural  transition  @  54  cm 
Delay  NOT  observed 


•  In  both  cases,  T*  and  Te  are  above  the  critical  960K  for  C02 

•  Why  is  delay  not  observed  for  injection  in  the  case  on  the  right? 
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Ongoing  Challenges  GAL$t 


Possible  Explanation 

*  N-factor  at  transition  is  similar  (5. 0-5.6  for  a  noisy  tunnel)  in  both  cases 

*  N-factor  at  the  injector  section  location  (13.3  cm  from  the  tip)  is  therefore  significantly 
higher  for  the  case  with  earlier  natural  transition  (right  hand  plot  on  previous  slide) 

*  To  suppress  the  2nd  mode,  mixing  must  be  achieved  at  a  relatively  low  (but  not  precisely 
known)  N-factor 

Variables  to  optimize  for  attaining  delay  by  injecting  C02 

*  N-factor  at  injection  location 

•  (T*  or  Te  of  boundary  layer  base  flow)  /  (Tvib  of  C02) 

*  (mixing  distance  for  C02  with  boundary  layer  base  flow)  /  (cone  length) 

•  (C02  mass  flow  rate)  /  (boundary  layer  base  mass  flow) 

Collaboration  with  Alexander  Fedorov  -  Moscow  Inst  of  Physics  and  Tech 

*  Determine  how  low  the  N-factor  at  injection  must  be,  and  where  this  physically 
occurs  on  the  cone 

*  Redesign  of  injector  section  to  move  it  closer  to  the  tip,  achieving  injection  before 
the  2nd  mode  acoustic  waves  appear 
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Conclusions 


GALM 


'  yitt^viw  ■  etrtiOP  I  MmvIA  * 


At  10  MJ/kg  enthalpy,  demonstrated  delay  versus  Argon  injection  and  also  versus  a  smooth 
injector 

—  C02  does  make  a  difference! 

Selected  a  new  condition  at  about  8-9  MJ/kg  for  further  study 

—  Meant  to  show  a  greater  effect  because  natural  transition  occurs  near  the  middle  of  the 
cone 

—  However,  C02  injection  did  NOT  seem  to  delay  transition  at  this  condition 
Designed  and  installed  a  half-porous,  half-smooth  porous  injector 

—  Provides  a  non-injection  “control”  with  every  injection  experiment 

Collaboration  with  Alexander  Federov  to  for  theoretical/computational  input  into  injector 
design  and  placement 
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Questions  ? 
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Back  up 
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Preliminary  results  from  resonantly  enhanced 
fjeldfocusedschlieren system  (REFFSS) 


GAL&t 


10/13/2011 
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Preliminary  results  from  resonantly  enhanced  field 
focused  schlieren  system  (REFFSS) 


10/13/2011 
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CALM 


Preliminary  results  from  resonantly  enhanced  field 
focused  schlieren  system  (REFFSS) 


10/13/2011 
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Preliminary  results  from  REFFSS 


RMS  Response  vs.  Time  -  Shot  2644 


10/13/2011 
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Preliminary  results  from  REFFSS  GALutt 


RMS  Response  vs.  Time  -  Shot  2647 


°iloo  1600  1700  1800  1900  2000  2100 

Time  (ps) 

h_R  =  6  MJ/kg  ;  p_R  =  30  MPa 


10/13/2011 
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Heated  Carbon  Dioxide  CALffit 


•  Baseline  condition  similar  to  shot  2541 

—  Test  gas  is  air 
—  Free-stream  Mach  is  5.3 
—  Isothermal  wall  at  293  K 

•  Pre-heated  C02 

—  Momentum  of  injection  matched  with 
13.5  g/s  of  cold  carbon  dioxide 

•  Increase  in  heating  results  in 
decreased  amplification 

—  Reduction  in  amplification  more 
efficient  near  1000  K 


Transpiration  interval 

25  \ 


Smooth  Cone 
T=  293  K 
T  =  1000  K 
T  =  1500  K 
T  =  2000  K 
T  =  2500  K 


\ 

Estimated  N  factor  oltranslion 


0.2  0.4  0.6  0.8  1 

s  (m) 


20  - 


15  - 


10  - 


5  - 


Stagnation  Conditions 

Free-stream  Conditions 

Pressure  (MPa) 

50.92 

Density  (kg/m3) 

0.05572 

Temperature  (K) 

5968.5 

Temperature  (K) 

1369.4 

Entalpy  (MJ/kg) 

9.51 

Velocity  (m/s) 

3957.9 

N  factor  versus  distance 
along  cone  surface 
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Computational  Model  CAL^ft 


•  Computations  done  using  STABL  software  suite 
—  Mean  flow 

•  2nd  order  accurate  fluxes 

•  Modified  Steger-Warming 

•  1st  order  Implicit  DPLR  method  for  time  integration 

•  Finite  rate  chemistry  and  T-V  energy  exchange 
—  Disturbances 

•  STABL  PSE-chem  solves  the  parabolized  stability  equations 

•  PSE  predict  amplification  of  disturbances 

•  Finite  rate  chemistry  and  T-V  energy  exchange 

•  Semi-empirical  eN  method  used  for  determining  transition  location 
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Stanton  number  (St)  [-) 


Transition  Determination 
Uncertainty 


GALffl 


Plot  of  St  vs  Rex  for  T5-2589;  P0  =  56.3  MPa,  h0  =  10.37  MJ/kg,  T0  =  6251.7  K 

-2 
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Transition  Reynolds  Number 


Porous  Injector  Results: 
Intermittency  Method 


calM 


Re  vs.  injection  @  ~10  MJ/kg,  ~55  MPa 


'  C02 
a  Argon 
Smooth 
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Intermittency  (7) 


Data  Reduction:  Intermittency  Method 


Shot  2600  Intermittency  (Average  by  Gauge  Position) 


Turbulent  intermittency 


0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1 


Shot  2600  F(y)  (Average  by  Gauge  Position) 


Transition  location 
determined  from  intersection 
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Intermittency  (7) 


Data  Reduction:  Intermittency  Method 


Shot  2596  Intermittency  (Average  by  Gauge  Position) 

0.25 

0.2 

0.15 

0.1 

0.05 

0 

0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1 

Shot  2596  F(y)  (Average  by  Gauge  Position) 


Turbulent  intermittency 


Transition  location 
determined  from  intersection 
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Stanton  number  (St)  [-] 


Porous  Injector  Results  (10  MJ/kg)  GAL&it 


10-micron  Porous  Injector  (Ar  injection  at  11.6  grams/sec) 


Plot  of  St  vs  Rex  for  T5-2600:  P0  =  54.7  MPa,  h0  =  9.88  MJ/kg,  T0  =  6044  K 


°  experimental 

laminar 

Van  Driest  II 

White  &  Christoph 

Reynolds  number  (Re^  [-] 

Transitional  flow 
Retr  =  2.88x106 


Heatflux  distribution  -  Shot  #2600 


Heat  Flux  (MW/m2) 
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Stanton  number  (St)  [-] 


Porous  Injector  Results  (10  MJ/kg) 


GAL&t 


10-micron  Porous  Injector  (no  injection) 


Plot  of  St  vs  Rex  for  T5-2598;  P0  =  55.3  MPa,  h0  =  10.26  MJ/kg,  T0  =  6204  K  Heatflux  distribution  -  Shot  #2598 


Initially  laminar  flow 


0.5  1  1.5  2  2.5  3  3.5  4  4.5  5 

Heat  Flux  (MVWm2) 


Retr  =  4.12  x  106 
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Stanton  number  (St)  [-] 


Porous  Injector  Results  (10  MJ/kg) 


GAL&t 


10-micron  Porous  Injector  (C02  injection  at  11.6  grams/sec) 


Plot  of  St  vs  Re^  for  T5-2590;  P0  =  55.9  MPa,  h0  =  9.92  MJ/kg,  T0  =  6064.3  K 


°  experimental 

laminar 

Van  Driest  II 

White  &  Christoph 

Reynolds  number  (Re^  [-] 

Completely  laminar  flow 


Retr  >=  5.22x10® 


Heatflux  distribution  -  Shot  #2590 
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Porous  Injector  Results  (10  MJ/kg)  GAL$tt 
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3  - 
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A  C02  (Porous  Injector) 

#  Smooth 

♦  Argon  (Porous  Injector) 


\ 


\  I 


1  \ 
\ 


J_ I_ I_ I_ I_ I_ I_ L 


0  2  4  6  8  10  12  14  16  18 

Injection  Mass  Flow  Rate(g/s) 

Summary  of  results  (average 
heat  transfer  method) 
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Transition  Reynolds  Number 


Porous  Injector  Results  (10  MJ/kg)  GALfflt 


Re  vs.  injection  @  ~10  MJ/kg,  ~55  MPa 


0  C02 
■  Argon 
A  Smooth 


mass  flux  (g/s) 


Two  porous  control  shots 
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Transition  Reynolds  Number 


Porous  Injector  Results  (10  MJ/kg)  GALfflt 


Re  vs.  injection  @  ~10  MJ/kg,  ~55  MPa 


0  C02 
■  Argon 
A  Smooth 


mass  flux  (g/s) 


Five  smooth  control  shots 
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Stanton  number  (St)  [-] 


Porous  Injector  Results  (6  MJ/kg) 


CALtit 


Plot  of  St  vs  Re^  for  T5-2582:  P0  =  29.6  MPa.  h0  =  S.69  MJ^kg.  T0  =  4129.1  K 


Porous  tip  -  no 
injection 
vacuum  plenum 


Plot  or  St  VS  Rex  for  15-2505;  R0  =  30.4  MPa,  h0  =  5.92  MJ/kg,  T0  =  4242.6  K 


Injection  of  C02  and  argon 
destabilize  boundary  layer 
and  transition  occurs  earlier 
at  lower  enthalpy/temp. 


106  107 
Reynolds  number  (ReJ  [-] 


Ar  injection  4  g/s 


C02  injection  4  g/s 
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Qualitative  visualization 


GALffl 
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Before  shot 


Porous  Injector  with  cone 
downstream  and  80  psi  injection 


Schlieren  seems  to  capture  the  C02  injection  clearly 
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Stanton  number  (St)  [-] 


Porous  Injector  Results 


GAL&t 


Heatflux  distribution  -  Shot  #2573 


Plot  of  St  vs  ReK  for  T5-2573;  P„  =  53  MPa,  h„  =  10.13  MJ/kg,  Tfl  =  6139.9  K 


o  experimental 

- laminar 

Van  Driest  II 
White  &  Christoph 


Reynolds  number  (Re  )  [-] 


C02  injection 
at  18.5  g/s 


Immediately 
turbulent  flow 
(streak  due  to 
injector  flaw) 


C02  injection 
at  26.6  g/s 
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I  i  I  T 
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Heat  Flux  (MW/m2) 

Heatflux  distribution  -  Shot  #2572 
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.5  5 
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Heat  Flux  (MW/nrT) 
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CFD  test  conditions  CAL&t 


Gas  Composition  (by 
mass  fraction) 

N2 

0.7345 

02 

0.1844 

NO 

0.0654 

N 

0.0 

0 

0.0157 
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AFRifi  Why  vibration  relaxation  and 
dissociation  damp  acoustic  wav 


Theory  known  for  decades  (Lighthill, 
1956,  Herzfeld  and  Litovitz,  1959, 
Clarke  and  McChesney,  1964, 
Vincentti  and  Kruger  1967) 


Following  Clarke  and  McChesney,  the 
linearization  of  perturbations  of 
the  N-S  equations  leads  to  damping 
curve  as  shown  -  Maximum  damping 
occurs  when  a)x=af/ae 


Relaxation  processes  such  as 
molecular  vibration  and  dissociation 
cause  damping  of  acoustic  waves 
through  phase  lag  between  pressure 
and  density 


oL 

10“' 
Equil.  -e 


10U 
(ojr„/a) 


10 

->  frozen 
From  Clarke  and  McChesney 
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injector  models 


Model  No. 

Hole  diameter, 
d  (mm) 

No.  Rows 

Injection 
angle,  a 
(deg) 

1 

0.51 

2 

12 

2 

0.76 

1 

11 

3 

0.76 

1 

12 

4 

1.02 

1 

11 

Model  No.  1 


Injector  Variants 


Model  No.  2 


Four  injectors  designed  and  built 
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Back-up  charts 


GALffl 
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Vibrational  temperatures  CAL&tt 


species 

A/ttj  o'©7’ot  h ^  qv  ©ti  qEq  gei  £e 

g/mol  K  J/mol  K  J/mol 

n2 

28.016  5.79  OO  I  3353.2  1  3  6.015  x  105 

6  7.136  x  10s 

1  7.342  x  105 

02 

32.000  4.16  OO  i  2239.0  3  2  9.225  x  101 

1  1.579  x  105 

3  4.320  x  10s 

3  5.960  x  10s 

Ar 

39.944  -  OO  -  -  I  5  1.115  x  106 

3  1.122  x  106 

N 

14.008  -  4.713  x  10s  -  -  4  6  2.301  x  lO* 

4  2.308  x  106 

6  3.452  x  10s 

12  9.971  x  10s 

O 

16.000  -  2.468  x  10D  -  -  5  3  1.903  x  103 

1  2.717  x  103 

5  1.899  x  10£ 

1  4.044  x  105 

5  8.829  x  105 

NO 

30.008  2.45  8.990  x  10*  1  2699.2  4  2  5.262  x  105 

4  5.496  x  10s 

C 

12.011  -  7.116  \  10c  -  -  I  3  1.962  x  102 

5  5.204  x  103 

5  1.219  x  105 

1  2.590  x  105 

5  4.036  x  10s 

co2 

44.011  1.13  -3.933  x  10^  2  960.1  1 

1  1992.5 

1  3380.2 

CO 

28.011  2.78  —1.139  x  lO11  1  3082.0  1  6  5.824  x  10b 

3  6.687  x  10s 

6  7.453  x  10s 

2  7.785  x  10s 

H2 

2.016  175.09  0.0  1  G  100.0  1  I  1.097  x  1QP 

2  1.196  x  106 

3  1.352  x  106 

H 

1.008  -  2.110  x  I0b  -  -  I  3  9.839  x  10* 

4  1.166  x  106 

5  1.230  x  106 

Oxygen  Carbon  Oxygen 

^  ^  f^\  Symmetric 

_ J  J  modes  (uOO) 


Bending 
modes  (OvO) 


(  A  _ g _ g  f  Asymmetric 

modes  (00u) 
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Dissociation  temperatures  CALffit 


W 


Air,  p  =  D.Olkg/'m3 


CO2,  p  =  0.01  kg/m’3 


N2,  p  —  0.01kg/m3  N2,  p  —  O.lOkg/m3 

St=q/rhoe*ue*(h0-0.5*ueA2*(1-r)-CpTw) 
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Detailed  condition  data 
(old  shots) 


calM 


Tkhle  A,  Li  Summary  of  freestream  conditions  far  all  shots. 


Shot 

Pc 

[MPa] 

To 

[K] 

ho 

[MJ/fcj] 

Poo 

[^Fa| 

[K] 

Poo 

[toj/m3] 

tfoo 

[m/a] 

l-l 

2131 

52.2 

7375 

10.75 

21.8 

1487 

4. 92X1C1-2 

4286 

5.6 

2332 

51.2 

7095 

10.90 

22.8 

1560 

4. 06xlO"2 

4211 

5.4 

2333 

50.0 

4434 

0.43 

40.9 

2214 

8.69x1 0"2 

324S 

4.3 

2334 

50.1 

5380 

10.89 

32.8 

1980 

6.29x10-2 

3709 

4.7 

2335 

37.6 

3459 

5,58 

27.7 

1426 

9.86xl0"2 

2647 

4.6 

2335 

37.2 

3606 

6,12 

27.8 

1532 

9. 12x1 0"2 

2741 

4.6 

2337 

51.2 

7379 

11.37 

23.2 

1662 

4.74X10"2 

4338 

5.4 

2433 

39.1 

5305 

3.34 

17.9 

1250 

5.00  xl0“2 

3712 

5.4 

2434 

42.5 

5538 

8.87 

20.3 

1360 

5.17x10-2 

3816 

5.3 

2435 

48.9 

5843 

9.53 

245 

1522 

5.60x10"  2 

3942 

5.2 

2435 

49.5 

6069 

10.07 

25.7 

1650 

5.39x1 0"2 

4040 

5.1 

2437 

46.8 

6567 

11.33 

27.3 

2178 

4.37x10"  2 

4363 

4.8 

2433 

48.9 

6069 

10.08 

25.5 

1660 

S3£xl0-2 

4040 

5.1 

2439 

55.5 

5715 

9.15 

26.6 

1423 

e^ixicr2 

3873 

5.3 

2440 

53.8 

5216 

7.99 

25.  L 

1182 

739xltr2 

3643 

5.4 

2441 

52.6 

46:89 

7.33 

27.7 

1285 

8. 05x1 0"2 

3407 

5.1 

2442 

53.3 

4250 

7.18 

33.0 

1479 

8.04xlO"2 

3201 

4.8 

2443 

53.9 

3702 

4.38 

36.3 

1443 

l.lSxlCr1 

2938 

4.7 

2444 

53.0 

3616 

5.99 

40.3 

1560 

l^lxlCr1 

274H 

4.5 

2445 

53.3 

3782 

6.61 

40.4 

1674 

mxlCr1 

2834 

4.5 

2446 

55.8 

4742 

6,85 

22.4 

931 

8.36x1 0"2 

3384 

5.7 

2447 

50.5 

4357 

6.06 

18.9 

783 

8.43x10"  2 

3193 

5.8 

2443 

48.4 

4358 

6.07 

13.0 

783 

8.01x10"  2 

3103 

5.8 

2449 

45.7 

4286 

5,94 

17.1 

76£ 

7.34x10-2 

3164 

5.8 
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Stanton  number  (St)  [-] 


h***  '  ■  etviiOP  i  mmvia 


Preliminary  Results  with  Porous  Injector^y^/#/ 


Control  Experiment  - 
Smooth  Surface 


Porous  injector  design 
looks  very  promising 


Ptiti  of  St  us  Fte_  lui  T5  2510:  P.  =  5t.fi  MPa.  K  =  9.47  MJukq,  T*  =  5052.fi  K 
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C02  Injection  System  GAL&tt 


INJECTOR 
MODEL 
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Experimental  model  GALfflt 


Injector  section  installed  in  ~1m  long,  5-degree  half  angle  cone 
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Stanton  number  (St)  [-] 


Preliminary  Results  (10  MJ/kg) 


GAL&t 


Control  Experiment  -  Smooth  Surface 


Plot  of  St  vs  Rex  for  T5-2609;  P0  =  55.7  MPa,  h0  =  9.92  MJ/kg,  T0  =  6063  K 

-2 


10 


Heatflux  distribution  -  Shot  #2609 


10 


Reynolds  number  (Re^)  [-] 

■0.3  ■ _ i _ _ i _ i _ i _ i _ i _ " 

-0.9  -0.8  -0.7  -0.6  -0.5  -0.4  -0.3 

Cone  surface  coordinate  (m) 

Initially  laminar  flow 

0.5  1  1.5  2  2.5  3  3.5  4  4.5  Z 

Heat  Flux  (MW/nT) 


Retr  =  4.36  x  106 
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Four  rows  of  orifices 


GALffl 


i  «H*cgrvfft  >  StrtiC*  i  Mm-rtft  ■■ 


Injector  with  four  rows  of  orifices, 
installed  in  T5  test  section 
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Stanton  number  (St)  [-] 


Preliminary  Results 


CALM 


Four-Row  Injector  (C02  injection  at  26.0  grams/sec) 


Healflux  distribution  -  Shot  #2520 


Cone  surface-  coordinate  (m) 


Immediate  Transition 


Heat  Flux  (MW/m2} 
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Stanton  number  (St)  [-] 


Preliminary  Results 


CALM 


Four-Row  Injector  (no  injection) 


Heatflux  distribution  -  Shot  #2520 


-D.0  HO.B  -D.7  -0.0'  -0.5  -0.4  -0.3 


Cone  surface  coordinate  (m) 


Immediate  Transition 


Heat  Flu*  (MWm2) 
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1 


One  row  injector 


L4JMB  I  CH*Cqrvfft  .  i  MmvIA 


•-  I  ^ 


■ 


One  row  of  orifices 


Distribution  A:  Approved  for  public  release;  distribution  unlimited 


62 


Stanton  number  (St)  [-] 


Preliminary  Results 


CALM 


One-Row  Injector  (no  injection) 


Heatflux  distribution  -  Eliot  #2536 


Transition  (though  not 
immediately  full  turbulence) 


Heat  Flux  (MVWm2) 
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N  Factor 


T  profile  and  N  factor  for  high  and  low  enthalpy 

conditions 


"cmM 


Shot  2562  •  No  i  ejection 
P0  -  29.6  MPa,  hc  *  5,69  MjMg 


Shot  2569  *  No  injection 


for  Po=30MPa, 
ho=5.7  MJ/kg 
(shot  2582) 


Approximate  transition 
location  -  42  cm 


0.2 


0,4  0*6 

s(m) 


Mach  Number 

4  6 


0,3 


10 


T :  solid  line 
M  :  dash  dot  line 


Mach  Number 
4  6 


10 
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Transition  Reynolds  Number 


N2  results 


Gas  inj  into  N2  @  ~10  MJ/kg,  “50  MPa 


5.00E+06 

4.50E-iO6 

4.00E+06 

3.50E-iO6 

3.00E+06 

2.50E4O6 

2.00E+06 

1.50E-iO6 

1.00E+06 


mass  flux  (g/s) 


GAL&t 


•  C02 
■  Argon 
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Porous  Injector  Results  (10  MJ/kg)  GALfflt 


x  106 


Summary  of  results 
(intermittency  method) 
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Transition  Reynolds  Number 


Porous  Injector  Results  1/2  CAL$t 


Re  vs.  injection  <§>  ~10  MJ/kg,  ~55  MPa 


C02 
■  Argon 
Smooth 


mass  flux  (g/s) 
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